Interstellar H2 toward HD 37903 
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Abstract 

We present an analysis of interstellar H2 toward HD 
37903, which is a hot, B 1.5 V star located in the NGC 
2023 reflection nebula. Meyer et al. (2010) have used a 
rich spectrum of vibrationally excited H2 observed by the 
HST to calculate a model of the interstellar cloud toward 
HD 37903. We extend Mayer's analysis by including the 
y"=0 vibrational level observed by the FUSE satellite. 

The Toi temperature should not be interpreted as a rota- 
tional temperature, but rather as a temperature of thermal 
equilibrium between the ortho and para H2. The ortho to 
para H2 ratio is lower for the lowest rotational levels than 
for the higher levels populated by fluorescence. The PDR 
model of the cloud located in front of HD 37903 points 
to a gas temperature T|;„= 110-377 K, hydrogen density 
n//= 1874-544 cm^-' and the star-cloud distance of 0.45 
pc. 

Key words: ISM: clouds — ISM: molecules — ultravi- 
olet: ISM 



1 Introduction 



A rich spectrum of vibrationally excited H? in the di 



rection to HD 37903 was first described by iMeyer et al. 



(|200lh . They have observed over 500 interstellar H2 ab- 
sorption lines from excited vibrational levels v"=l-14 and 
rotational levels up to T'=13. These lines were detected 
in a Hubble Space Telescope (HST) spectrum made with 
the Space Telescope Imaging Spectrograph (STIS). A Far 
Ultraviolet Spectroscopic Explorer (FUSE) spectrum was 
made after Mayer's publication allowing to access the 
y"=0 vibrational level of the ground electronic stat e. The 
FUSE spectrum was used by lRachford et a/.l(l2009h to de- 
termine the Top = \1QKI ln(9A^o/M)-68 ± 7 K gas "ki- 
netic" temperature and the hydrogen molecular fraction 
f(H2)=2N(H2)/(2N(H2)-FN(HI))=0.53 ± 0.09 in the direc- 
tion towards HD 37903. 

The star HD 37903 was also observed by the Berkeley 



Extreme and Far-Ultraviolet Spectrometer (BEES) on- 
board the Orbiting and Retrievable Far and Extreme Ul- 
traviolet Spectrometer ( OREEUS) telescop e. The spectral 



resolution was R=3000. Lee et al. (20021) have used this 



spectra to obtain H2 column densities on y"=0 and J"=0- 
5 rotational levels. Their column densities agree within 
the order of magnitude to the column densities derived 
in this paper from the FUSE spectr a. The physical pa- 
rameters derived by Lee et al. ( 200 2b are: Top=63 ± 5 K, 
f(H2)^0.496 ± 0.017 and the cloud density n=5600 cm'l 

The H2 molecule exists in two forms: ortho (odd J") 
and para (even J") H2. It is caused by the spins of the 
hydrogen nuclei which can point "in the same" direction 
(ortho H2 - triplet state, total nuclei spin 1=1) or in op- 
posite directions (para H2 - singlet state, 1=0). The ratio 
ortho/para H2 is therefore 3:1 at standard temperature and 
pressure. Conversion between this two spin isomers can 
take place in gas phase, or on the surface of gas grains 
(ILe Bourlotil2000r) . The ortho-para conversion in the gas 
phase is caused by the exchange of proton in collisions 
with H, H+ and H+. 

The fluorescence cascade of H2 that leads to popula- 
tion of excited ro-vibrational states of H2 as well as to 
the emission of infrared phot ons from quadrupo l e tran - 
sitions has been described by iBlack & Dalgamo 
The first ultraviolet detection o f vibrationally excited in- 
terstellar H2 was performed by iFederman et al.\ (Il995h in 
the HST spectrum of ^ Ophiuchi. 



2 Column densities 

We have used both HST STIS and FUSE spectrum to ob- 
tain column densities on H2 ro-vibrational levels. The 
HST STIS spectrum o59s04010 was averaged because it 
consist of two subexposures. The HST spectrum is lo- 
cated between 1160 and 1357 A. We have also used the 
FUSE observation PI 160601. We have analyzed only 
spectrum from detectors IB LiF and lA LIE. These two 
spectra had the best quality. The FUSE spectra originat- 




Figure 1 : A fragment of FUSE spectrum with a fit (green line) of H2 absorption lines. Tlie red line represents 
continuum. 
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Figure 2: Fragment of HST STIS spectrum (gray line with dots) fitted with 268 H2 absorption lines (green line). The 
figure presents only a 5 A fragment of the spectrum, while the presented fit was done to the whole 200 A long STIS 
spectrum (1160-1357 A) and included 7449 H2 lines. 
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Figure 3: Occupation of H2 X ro-vibrational levels towards HD 37903. Note the different y-axis scale for v"=0 (left) 
and for higher y" levels (right). 



ing from the same detector were shifted and coadded us- 
ing the IRAF tasks pojfsets and specalign. A part of the 
FUSE spectrum is shown on fig. [T] The whole FUSE 
spectrum ranges from 987 A to 1188 A. 

The Ar I 1048 A Une that lies in the wing of R(0) and 
R(l) lines was cut out from the spectrum. The wings of 
H2 lines were calculated up to 60 A from the line center. 
The absorption lines in the FUSE spectrum were modeled 
with the Voigt function, while a Gauss function was used 
for modeling the H2 lines in the STIS spectrum. 

We have used a Gaussian poin t spread function (P SF) 
with FWHM equal to 15 km/s djensen et al\ boioh for 
modeling the FUSE spectra. The column densities at 
J"=0, 1 and 2 (v"=0) were derived from transitions be- 
tween B and X electronic levels. The vibrational transi- 
tions used are (0,0), (2,0), (3,0) and (4,0), where the first 
digit is the v' on the upper electronic level B and the sec- 
ond digit v" = is the vibrational level of the ground 
electronic state X. The FUSE spectrum at the (1,0) vibra- 
tional transition was to noisy to perform a good fit. The 
H2 li ne positions and osc illator strengths were adopted 
from lAbgrall et al. (1994). The total transition probabil- 
ities that include the transi tions to continuum ( dissociat- 
ing) states were taken from lAbgrall et al. ( 2000 ). 

The column densities of rotational levels J"=0-9 from 
the vibrational v"-0 level as well as the column density 
of the v"=l J"=0 level were derived from the FUSE spec- 
trum. Other ro-vibrational levels of the X ground elec- 
tronic state were derived from the o59s04010 HST STIS 
spectrum. 

The column densities were derived with the profile fit- 
ting technique. At each point of the simulated spectrum 
the optical depth of many spectral lines has been summed. 
Such procedure was necessary to calculate the profile of 
blended lines which were ubiquitous in the spectra. The 
cloud velocity, the doppler broadening parameter and col- 
umn densities on all observed levels were free parameters, 
that were fitted to the observed spectra. 

The STIS spectrum was fitted with H2 absorption lines 
from all vibrational v"=0-14 levels, and rotational levels 
J"=0-13. H2 lines that were blended with atomic lines 
were excluded from the fitting procedure. Total 7449 H2 
lines were included in the simulated STIS spectrum. The 
whole STIS spectrum (200 A long) was fitted at once with 
all 7449 H2 lines, because of large number of blended 
lines. A fragment of the HST STIS spectrum with lines 
of vibrationally excited molecular hydrogen is presented 
on Fig. [2] The observed column densities are presented 
in Table [U and on Fig. [3] The errors of the column den- 
sity are about 20% for J"=0-7 rotational levels, and up to 
~40% for higher rotational levels. 
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Figure 4: Theoretical N(J^1)/N(J=0) ratio on the v"^0 
vibrational level from eq. |6]is shown as a dashed blue line. 
The light blue dashed lines represent errors introduced by 
the uncertainty of the ^orthof^para H2 ratio 0.64 ± 0.21. 
The observed N(1)/N(0) ratio together with its uncertainty 
is plotted as a straight line. 

3 Results 

Total observed column density of H2 on all vibrational 
and rotational levels equals to N(H2)=(9.6 ± 1.6) • 10^" 
cm~^. The observed H2 column densit y is slightly 
higher than N(H2)=4- 10^° cm"^ predicted by lMeyer et al. 
(|200lh . Co lumn density of neutral hydrogen N(HI)= 1 .48 ■ 
10 cm" doiplas & Savagelll994h . The hydrogen molec- 
ular fraction (assuming 10% error of N(HI)) equals to 
f(H2)-2N(H2)/(2N(H2)-FN(HI))=0.56±0.04. 

The Top temperature, calculated from the lowest en- 
ergy levels of ortho and para H2 (v" = 0) was derived 
from the equation: 



exp - 



A^(0) gmgAO) 



IcTop 



(1) 



where gj is the spin degeneracy factor g[{I) = 21+1. 
The gi{l) = 3 is the statistical weight for ortho H2, and 
giiO) = 1 for the para-H2 spin izomer. The. gj{J) = 27-1-1 
is the statistical weight for the rotational level J. The Tqp 
is the temperature of thermal equilibrium between the or- 
tho and para spin izomers. The derived Top=67±8 K is 
simila r to the temper atures determine d by lRachford et al. 
(l2009h (68 K) and bv'Le egfa/.! t00± (63 K). The photon 
dominated regions (PDR) models of interstellar clouds 
shows that the Tgp t emperature is correl ated with the gas 
kinetic temperature JUe Petit et al. 200 ^. 

However, if we want to calculate the rotational temper- 
ature across the ortho-para divide we have to take into 
account the ^ortiwl^para H2 ratio. The population of the 
ro-vibrational levels of H2 depends not only from the tem- 



Table 1: Column densities of the H2 ro- vibrational levels towards HP 37903 [cm 
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Figure 5: Rotational temperatures for para and ortho H2 
toward HD 37903 as a function of the vibrational number 
y". 



Table 2: Rotational temperatures for para and ortho H2 
toward HD 37903. 
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perature and radiation field, but also from the total ortho 
to para ratio. In order to include the ortho to para ratio 
into the Boltzmann equation we will write separate Boltz- 
mann distributions for the ortho and para spin izomers. 
The Boltzmann distribution for the ortho H2 can be writ- 
ten as: 



E{Jo") 
kT 



(2) 



where N ortho is the total amount of ortho H2. The partition 
function Z is: 



(3) 



Let us assume, that para H2 has a Boltzmann distribu- 
tion with the same rotational temperature: 



NpiJp") = 



N para 



g/(0)g7(-/p")exp 



EjJp") 
kT 



(4) 



where Npara is the total amount of para H2, and the parti- 
tion function for para H2: 



2] gmgAJp")expl- 

Jn" even ^ 



EjJp") 
kT 



(5) 



By dividing the equations |2]by|4] we obtain: 

NoiJol N ortho gl{l)gj{Jo") Zp{T) 



NpiJp") 



Npara gliO)gjiJp")ZoiT) 
E(Jo") - EiJp") 



■ exp 



kT 



(6) 



We have tried to solve the above equation numerically for 
the Jo"=l and Jp"=0 states (y"=0) to obtain the Tqi rota- 
tional temperature. The observed Northo /Npara - 0.64 ± 
0.21. Figure IHpresents the right side of eq. [6]and the ratio 
of observed column densities (left side of eq. O together 
with maximal errors. The rotational temperature Tqi can 
take values between and 500 K. 



The equation |6] can be written in another way by mak- 
ing two assumptions that are not valid in the interstellar 
medium: 



1. r>240K = 

2- A'^ ortho I N para 

tions. 



Zp{T)IZ„{T) « 1/3 

= 3 as in standard laboratory condi- 



Under this assumptions the eq. [6] for rotational tempera- 
ture reduces to: 

■exp — |. (7) 



NpiJp") gmgAJp") 



kT 



and Tq\ is the same as the temperature Tqp of ortho-para 
thermal equilibrium. 

The rotational temperatures T02 and T13 involve ro- 
tational levels of the same spin isomers of H2, and do 
not depend on the total ^orthof^para H2 ratio and on the 
partition function. For the cloud towards HD 37903 the 
T02 = 132 ± 14 K and Tn = 172 ± 9 K. However, the T12 
temperature calculated from eq. [6] is 1 33^^25 while the 
use of eq. |7]gives a incorrect temperature of 263 + 89 K. 

A good illustration that the eq. [6] must be used can be 
the y=l level. On the v - \ level the temperature T12 = 
1223 K (levels ortho and para) is close to the T02 = 921 
K and T13 - 1090 K temperatures calculated from the 
same spin isomers. The temperature calculated with the 
traditionally used equation |7] is negative T= -466 K ! 

The ortho (odd J") to para (even J") H2 ratio (here- 
after 0/P) (0//')y"=i_i4 = 1.35 ± 0.18 was calculated by 
summing ortho and para H2 on vibrational levels v"-\- 
14. It is a bit lo wer than 0/P=1.45 ± 0.08 given by 
Meyer et al\ (l200lh. We have also used the method pro- 
posed by Iwilgenbus etaU (l2000l) to calculate 0/P H2 
from individual ortho levels, specially 7o"=l and 3 on 
y"=0. For the y"=0 J"=l ortho state the 0/P ratio dif- 
fers significantly from the (C?/P)v"=i-i4 = 1.35 value 
and is equal to {0/P)j"=i = 0.63 ± 0.11. Also for the 
y"-0 J"=3 ortho state the observed 0/P is low - only 
{0/P)r=3 = 0.70 ± 0.07. For higher rota tional and vi- 
brational levels the Wilgenbus et al. ( 2000h method gives 
(0/P)/^, ratio close to 1.3 derived for excited vibrational 
H2 states. 

The rotational temperatures for vibrational levels 
y"=l-14 for para H2 were obtained from the linearized 
Boltzmann distribution: 



hi 



NpiJp") , ^^^^^ E{Jp") 
gi(0)gAJp ) kTy 



(8) 



Analogous equation was used for the distribution of the 
ortho H2 states. Fig. [3] shows left side of eq. [8] versus 
E/k. Levels fulfilling the Boltzmann distribution should 



Table 3: Comparison of cloud models toward HD 37903. 
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be placed on a straight line on this plot. On the y"=0 
vibrational level, which is populated partially by colli- 
sions the occupation of rotational levels does not follow 
a straight line. Therefore the rotational temperature was 
not calculated for the y"=0 vibrational level. For both spin 
isomers and for vibrational levels y"=l-14 the tempera- 
ture was calculated using the the linear regression method. 
The inverse of the line inclination {-l/Ty in eq. (8]) gives 
us the rotational temperature. The resulting temperatures 
are presented in table [2] and on figure [5] The rotational 
temperatures for the ortho isomer are usually higher than 
for the para H2. However the ortho H2 rotational temper- 
atures calculated form levels Jo=l to last but one ortho 
state are in conformity with the para H2 rotational temper- 
atures. So the assumption that the rotational temperatures 
Ty are the same for ortho and para H2 seems to be jus- 
tified. The rotational temperatures falls hnearly with the 
increasing vibrational level. 



4 Model 



We h ave used the Meudon PDR code (ILe Petit et al. 
'200^ to calculate a model of the interstellar medium in 
the direction of HD 37903. The interstellar reddening 
E(B-V)=0"? 35 and Rv=5 . 5 was chosen to match the val- 
ues used by iMeyer et all (120011) . The radiation source is 
the HD 37903 star (B 1.5 V spectral type). The model in- 
cludes 300 H2 ro-vibrational levels and uses exact compu- 
tation of radiative transfer in H2 spectral lines. The model 
is isobaric with thermal balance. 

All our models were two-side models with interstellar 
radiation field on the observer side equal to one 'Draine' 
unit. The interstellar radiation field on the star side has 
been varied from 1 to lO'* 'Draine' units. The best model 
was chosen among 6917 different models by minimising 
the sum: 



2] Wy\r{loi 

v"J" 



NobsiO,0) 



log 



Nmodeliv", J") 



N, 



model 



(0,0) 



(9) 



where the weights Wy'^r were chosen so, that the levels 
populated by collisions J"=0-6 (v"=0) have the same in- 



fluence on the final sum as the rest of the levels (v">0 and 
v"=0 J"=7-13), that are populated by fluorescence. 

We have varied the star-cloud distance, hydrogen den- 
sity and the interstellar radiation field on the star side in 
order to find a model that matches the observations. 

Table [3] presents of comparison of our best model and 



model s presented by iMeyer et al.\ (120011) and iLee et al. 
(|2002h . The star - cloud distance of 0.45 pc, that is 



responsible for the filling the fluorescence H2 levels is 
similar in our mo del and in the model presented by 
Meyer et al.\ (120011 ) (d=0.5 pc). Our cloud kinetic temper- 
ature (connected with the collisional-filled levels) equals 
to Tjiin-'ill K on the star side of the cloud and Tfa„=110 
K on the observer's side. The hydrogen density obtained 
in our best model changes from n//=544 cm"^ on the star 
side to 1874 cm"^ on the observer's side. The interstellar 
radiation field on both sides of our best cloud model is 
one 'Draine' unit. 



5 Conclusions 
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Column density of molecular hydrogen toward HD 37903 
on all observed H2 levels N(H2)^9.6 • 10^° cm^^ jhe 
hydrogen molecular fraction f(H2)=0.56. The ortho/para 
H2 ratio equals to 1.35 for the excited levels, but for the 
J"=l ortho state (v"=0) the ortho to para ratio is only 

(O/P)y=i=0.63. 

The Top is a temperature of thermal equilibrium be- 
tween the ortho and para spin isomers and is equal to 67 
K. The rotational temperatures To2=132 K and Ti3=172 
K. The formula for rotational temperatures calculated 
across the ortho-para divide (like T12) should include the 
(0/P)y H2 ratio. 

The best PDR model for the cloud toward HD 37903 
gives a gas kinetic temperature T<:,>,= 110-377 K, hydro- 
gen density n//=544-1874 cm"^ and star - cloud distance 
of 0.45 pc. 
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